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Abstract. Absorption spectra of small molecular aggregates (oligomers) arc 
considered. The dipole-dipole interaction between the monomers leads to shifts of the 
oligomer spectra with respect to the monomer absorption. The line-shapes of monomer 
as well as oligomer absorption depend strongly on the coupling to vibrational modes. 
Using a recently developed approach [Rodcn et. al, PRL 103, 058301] we investigate 
the length dependence of spectra of one-dimensional aggregates for various values of 
the interaction strength between the monomers. It is demonstrated, that the present 
approach is well suited to describe the occurrence of the J- and H-bands. 



1. Introduction 

Due to their unique optical and energy-transfer properties molecular aggregates have 
attracted researchers for decades [U EJ |3]. In recent years there has been a growing 
interest in one- dimensional molecular aggregates consisting of a small number N 
(ranging from only two to a few tens) of molecules [H EJ EJ El EJ EJ [10], where finite 
size effects play an important role. Such systems are e.g. dendrimers [8], stacks of 
organic dye molecules [5] or the light harvesting complexes in photosynthesis [HI (121 EE] . 
Usually in these aggregates there is negligible overlap of electronic wavefunctions 
between neighboring monomers. However, there is often a strong (transition) dipole- 
dipole interaction, that leads to eigenstates where an electronic excitation is coherently 
delocalized over many molecules of the aggregate, accompanied by often drastic changes 
in the aggregate absorption spectrum with respect to the monomer spectrum. Since 
these changes depend crucially on the conformation of the aggregate, commonly the first 
step in the investigation of such aggregates is optical spectroscopy. Important properties, 
like the strength of the dipole-dipole interaction can often already be deduced from these 
measurements (H [H HH [15]. For small aggregates, often referred to as oligomers, the 
finite size leads to pronounced effects in their optical properties. 

The interpretation of measured absorption spectra is complicated by the fact that 
the electronic excitation couples strongly to distinct vibrational modes of the molecules 
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and to the environment, leading to a broad asymmetric monomer absorption spectrum 
that often also shows a vibrational progression [U [TJ]. Because of this (frequency 
dependent) interaction with a quasi-continuum of vibrations, the theoretical treatment, 
even of quite small oligomers (of the order of five to ten monomers), becomes a 
challenge. Exact diagonalization [HI EJ EJ [T71 [18] is restricted to a small (one to five) 
number of vibrational modes. Furthermore, to compare with experiments, the resulting 
"stick-spectra" have to be convoluted with some line-shape function to account for 
the neglected vibrations. For large aggregates the "coherent exciton scattering" (CES) 
approximation, which works directly with the monomer absorption line-shape as basic 
ingredient [TH], has been shown to give very good agreement with experiment [T^l 1201 12T] . 
However, one expects problems for small (N < 10) oligomers [6]. 

To treat oligomers with strong coupling to vibrational modes, recently an efficient 
approach was put forward [22] which is based on a non-Markovian stochastic Schrodinger 
equation. In the present paper this approach will be used to investigate the absorption 
of linear one- dimensional oligomers. 

The paper is organized as follows: In section [2] the Hamiltonian of the aggregate is 
introduced followed by a brief review of our method to calculate the absorption spectrum 
in section [3j The calculated oligomer spectra are presented and discussed in section HJ 
We conclude in section [5] with a summary and an outlook. 

2. The model Hamiltonian 

In this work we will use a commonly applied model of a molecular aggregate including 
coupling of electronic excitation to vibrations [9j [231 121] • We consider an oligomer 
consisting of N monomers. The transition energy between the electronic ground 
and electronically excited state of monomer n is denoted by e n . We will investigate 
absorption from the state where all monomers are in their ground state. A state in 
which monomer n is excited and all other monomers are in their ground state is denoted 
by 1 ir n } . States in which the whole aggregate has more than one electronic excitation 
are not considered. We expand the total Hamiltonian 



of the interacting monomers and the vibrational environment with respect to the 
electronic "one-exciton" states | ir n ). The purely electronic Hamiltonian 

N 



contains the interaction V nm which describes exchange of excitation between monomer 
n and m. The environment of each monomer is taken to be a set of harmonic modes 
described by the environmental Hamiltonian 



H = H c \ + H int + H ( 



cnv 



(1) 




(2) 



N 




(3) 



71=1 j 
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Here a n j denotes the annihilation operator of mode j of monomer n with frequency oo n j. 
The coupling of a local electronic excitation to its vibrational environment is assumed 
to be linear and the corresponding interaction Hamiltonian is given by 

N 

H int = -^2Yl Kn i A + a nj)\^n)(^n\, (4) 
n=l j 

where the coupling constant K n j scales the coupling of the excitation on monomer n 
to the mode j with frequency u n j of the local vibrational modes. This interaction is 
conveniently described by the spectral density 

J n {u) = ^ \K n .j\ 2 5(uJ - L0 nj ) (5) 

i 

of monomer n. It turns out that the influence of the environment of monomer n is 
encoded in the bath correlation function [23], which at temperature T reads 

du J n (u>) ( cos(W) coth — — — — i sin(wr) J . (6) 
V 2k B T / 

In this work we restrict the discussion to the case T = in which the bath-correlation 

function a n {r) reduces to 

a n {r) = j du J n {uj)e- iUT = ^ M V^. (7) 

When the environment has no "memory", i.e. a n (r) oc 5(r) it is termed Markovian, 
otherwise it is non-Markovian. 



3. Method of calculation 



We use a recently developed method to describe the interaction of the excitonic system 
with a non-Markovian environment [22] which is based on ideas from a stochastic 
Schrodinger equation approach to open quantum system dynamics j25j [26] . Within 
this approach it is possible to exactly describe the dynamics governed by a Hamiltonian 
of the form of Eq. ([T]) by a stochastic Schrodinger equation in the Hilbert space of the 
electronic system alone. However, due to the appearance of functional derivatives in 
this stochastic Schrodinger equation, the general method is of limited practical use. To 
circumvent this problem in Ref. [22] these functional derivatives were approximated by 
operators D^ m '(t, s) which are independent of the stochastics and can be obtained from 
a set of coupled differential equations. 

Within this approximation the cross-section for absorption of light with frequency 
Q at zero-temperature is simply given by 

a{fl) = —QRe dte l ™(^(0)|^(t)) (8) 
"< c Jo 

with initial condition 

N 

Mo)> = £(£-A0k»>, ( 9 ) 

n=l 
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Figure 1. Sketch of the geometry considered for the case N = 5. The points mark 
the positions of the monomers and the arrows indicate the direction of the respective 
transition dipolcs. 



where the orientation of the monomers enters via the transition dipoles jl n and the 
polarization of the light S. The state \ijj{t)) is obtained from solving a Schrodinger 
equation in the small Hilbert space of the excitonic system 

dt\m) = - \H cl m)) + J2 Wm)(-Km\D( m) (t)m)) (io) 

m 

with initial condition fl9]). The operator D^ m \t) appearing on the right hand side of 
equation ffTUj) contains the interaction with the environment upon excitation of monomer 
m in an approximate way [22], and is given by 

£>M(t) = fdsa m {t-s)D^ m \t,s), (11) 
Jo 

where a m (t — s) is defined in Eq. (13) and the operator D^ m '(t, s) is obtained by solving 

d t D^ m \t,s)= [- l -H el ,D^\t,s)\ (12) 

+ J2[\^)(^\D {l) (t),D^\t,s) , 
i 

with initial condition D^ rn \t = s,s) = — |7r m )(7r m |. 

This method allows a numerically fast and simple treatment of an assembly of 
interacting monomers with a complex environment. 



4. Results 



Using the method described in the previous section we will now investigate the 
dependence of oligomer absorption spectra on the number N of monomers and 
the interaction strength between them. To this end we choose a one dimensional 
arrangement of the monomers and take the transition energies to be equal, e n = e. 
We take only coupling between nearest neighbors into account, which we assume to be 
equal for all neighboring monomers, i.e. V n ^ n+ i = V for all n. For simplicity we also 
take all transition dipole moments /2 n to be identical and parallel to the polarization of 
the light. This arrangement is sketched in Fig. [U Before including the complexity of 
the vibrations we will briefly discuss the purely electronic case, given by Eq. [2j For the 
arrangement described above the Hamiltonian Eq. [2] can be diagonalized analytically to 
obtain the eigenenergies Ej = e + 2Vcos(7T7'/ (N + 1)) where the sign and magnitude of 
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Figure 2. Absorption spectrum of a single monomer taken in the calculations. The 
inset shows the corresponding spectral density (unit 10 -2 eV/H). 



V depend crucially on the angle a between the transition dipoles and the axis of the 
oligomer. The oscillator strength Fj for absorption from the electronic ground state of 
the oligomer to an excited state with energy Ej is given by f27\ [T5| |2"8"] 

Fi = -ctg 2 — — r 13 

3 N + 1 & 2(iV+ 1) K J 

where the oscillator strength of a monomer is taken to be unity. From this one sees 
that the state with j — 1, which is located at the edge of the exciton band, carries 
nearly all the oscillator strength. For large N it carries roughly 81% of the oscillator 
strength; for N < 7 even more than 90%. Depending on the sign of V this state is either 
shifted to lower energies (for V < 0) or to higher energies (for V > 0) with respect to 
the electronic transition energy e. Note that with increasing number of monomers N 
this state is shifted further away from the monomer transition energy, approaching for 
iV — > oo the value E\ = 2V. Note also that the absorption spectra for V and — V 
are just mirror images of each other. This will change drastically when one includes 
coupling to vibrations. 

For the vibrational environment we consider a continuum of frequencies so that 
the spectral density J n {oS) becomes a smooth function. Exemplarily we will consider 
the monomer spectrum shown in Fig. [2j This spectrum is obtained from a spectral 
density (shown in the inset of Fig. [2]) which is composed of a sum of Lorentzians. The 
zero of energy is located at the transition energy e. Due to the strong interaction with 
the vibrations this spectrum is considerably broadened. It also shows a pronounced 
asymmetry (it is much steeper at the low energy side). This strong asymmetry is mainly 
due to coupling to high energy vibrations with vibrational energies around 0.18 eV, which 
is roughly the value of a C=C stretch mode in organic molecules. 

We will now take a closer look at the dependence of the oligomer spectra on the 
number N of monomers for various values of the interaction strength V. Each column 
of Fig. [3] shows, for a distinct value of the dipole-dipole interaction V, a sequence of 
absorption spectra with increasing length N of the oligomer from top to bottom. The 
values of V and N are indicated in the individual plots. On top of each column the 
monomer spectrum from Fig.[2]is shown again for comparison. In addition to the spectra 
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Figure 3. Absorption spectra of linear oligomers for different numbers of monomers 
N and various coupling strength V (in eV). The solid curves are calculated using the 
method described in section [3] for the spectral density of FiglJJ sticks are obtaind from 
the purely electronic Hamiltonian. 



calculated including all vibrational modes also "stick-spectra" resulting from the purely 
electronic theory are shown. 

The values chosen for the interaction strength V between neighboring monomers 
are V = ±0.15 eV (column 1 and 4) and V = ±0.05 eV (column 2 and 3). These values 
are of the order of magnitude as what is found in experiment [HE]. 

As explained above, the purely electronic stick spectra for opposite sign (but same 
magnitude) of V are mirror images with respect to each other. However, as can clearly 
be seen, this is no longer the case for the spectra including coupling to the vibrations. 
Although, in accordance with sum rules [T9l [29| 124"] . the mean of each absorption 
spectrum is centered at the mean of the corresponding purely electronic spectrum, the 
line-shape is completely different for negative and positive interaction V. 

We will now discuss the dependence of the absorption spectra on iV and V in 
more detail. For large negative interaction V = —0.15 eV (first column) the absorption 
line-shape considerably narrows upon increasing the number N of monomers of the 
oligomer. This exemplifies the appearance of the famous so-called J-band [2], a very 
narrow red-shifted peak found in the absorption spectra of large dye- aggregates. It 
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shows a strong asymmetry, with a steep slope at the low energy side and a long 
tail extending to higher energies. Due to this asymmetry of the line-shape the peak 
maximum appears shifted to lower energies with respect to the purely electronic "stick" . 
For weaker interaction V = —0.05 eV (second column) the narrowing can already be 
seen, but is less pronounced than for V = —0.15 eV. Furthermore, when increasing the 
number of monomers from 3 to 6 there is little change in the absorption spectrum for 
V = —0.05 eV, while for V = —0.15 eV there is still a noticeable narrowing. In contrast 
to the narrowing found for negative values of V, for the same absolute value \V\ but 
positive V, the spectra are much broader and show little resemblance with their negative 
counterparts. This is the typical case of the so-called H-band, which is blue-shifted with 
respect to the monomer absorption. 

Note that for the present values of V there is only minor change in the absorption 
spectrum for N > 6, which indicates that for larger N finite size effects will no longer 
be of great importance. 

5. Summary and Outlook 

In this work we have investigated the properties of one-dimensional linear oligomers 
using a recently developed approach based on a stochastic Schrodinger equation [22] 
to handle coupling of the electronic excitation to vibrations. Considering the length 
dependence of the oligomer spectra as well as the dependence on the dipole-dipole 
interaction between the monomers, we found that this approach, although approximate, 
is well suited to describe the basic features like the narrow J-band and the broad H-band. 
While in this work discussion was restricted to zero temperature, we plan to study the 
temperature dependence with the present method and compare to experiment. 
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